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ABSTRACT 

We analyze the unconventional top quark decay mode t H^b at the quan¬ 
tum level within the context of general Two-Higgs-Doublet models by includ¬ 
ing the full electroweak effects from the Yukawa couplings. The results are 
presented in the on-shell renormalization scheme with a physically well mo¬ 
tivated definition of tan 13. While the QGD corrections have been taken into 
account in the current experimental analyses of that decay, the electroweak 
effects have always been neglected. However, we hnd that they can be rather 
large and could dramatically alter the interpretation of the present data from 
the Tevatron collider. For instance, in large portions of the parameter space 
the electroweak effects prevent the Tevatron data from placing any bound at 
all to the charged Higgs mass for essentially any value of tan (3. 




With the discovery of the top quark at the Tevatron]^ the last matter building block of 
the Standard Model (SM) has been fully accounted for by experiment. Still, for an effective 
experimental underpinning of the fundamental mass generation mechanism in the SM one 
has to hnd the elementary Higgs scalar, which has intriguingly evaded all attempts up to 
now. Therefore, in spite of the great significance of the top quark discovery, the Higgs 
mechanism - the truly theoretical core of the SM - remains experimentally unconhrmed. 
On the other hand, the recent evidence for the possibility of neutrino oscillations |@] gives 
further support to the idea that the SM could be subsumed within a larger and more 
fundamental theory. The search for physics beyond the SM, therefore, has to continue 
with strong effort both at low and high energy. And, complementary to the low-energy 
experiments, the peculiar nature of the top quark - its large mass and its characteristic 
interactions with the scalar particles - may help decisively to unearth further vestiges of 
physics beyond the SM. 


The Two-Higgs-Doublet Model (2HDM)P plays a special role as the simplest exten¬ 
sion of the electroweak sector of the SM. In this class of models one enlarges the scalar 
sector of the SM by the introduction of another Higgs doublet, thus rendering the scalar 
sector in terms of 


After spontaneous symmetry breaking one is left with two CP-even (scalar) Higgs bosons 
/i°, a CP-odd (pseudoscalar) Higgs boson A° and a pair of charged Higgs bosons 
The parameters of these models consist of: i) the masses of the Higgs particles, M^o, 
Mfjo, Mjifi and Mh+ (with the convention M/jO < M^o), ii) the ratio of the two vacuum 
expectation values 


tan/3 = 



Us 


( 2 ) 


and the mixing angle a between the two CP-even states. Two types of such models have 
been of special interest 0 which avoid potentially dangerous tree-level Flavour Changing 
Neutral Currents: In Type I 2HDM only one of the Higgs doublets is coupled to the 
fermionic sector, whereas in Type H 2HDM each Higgs doubled {Hi, H 2 ) is coupled to 
the up-type fermions and down-type fermions respectively, the Yukawa couplings being 

g 1/2 Mw sin (3 ’ ^ g v^Mvk {sin/3, cos/3} 

Type H models do appear in specific extensions of the SM, such as the Minimal Super- 
symmetric Standard Model (MSSM) which is currently under intensive study both theo¬ 
retically and experimentally. 


In case that the charged Higgs boson is light enough, the top quark could decay via 
the non-standard channel t b. Based on this possibility the CDF collaboration at 
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the Tevatron has undertaken an experimental program which at the moment has been 
used to put limits on the parameter space of Type II models The bounds are obtained 
by searching for an excess of the cross-section a{pp —>■ ttX —>■ TVrX) with respect to 
a{pp —>• ttX —>• luiX) {I = e, fi). The absence of such an excess determines an upper 
bound on r(t —> b ^ Ur b) and a corresponding excluded region of the parameter 
space (tan/?, However, it has been shown that the one-loop quantum corrections 

to that decay width can be rather large. This applies not only to the conventional QCD 
one-loop corrections^ - the only ones used in Ref. |Q - but also to the QCD and elec- 
troweak corrections in the framework of the MSSM^, 0. Thus the CDF limits could 
be substantially modified by radiative corrections]^] and in some cases the bound even 
disappears. 

To our knowledge, and in spite of some existing approximate calculations Q, a fully- 
fledged account of the main electroweak corrections to r(t —> H^b) in the framework of 
general 2HDM’s is lacking in the literature. Thus, in this letter we address the complete 
computation of the one-loop electroweak contributions (EW) at leading order in both Aj 
and Af, in generic Type I and Type II 2HDM’s and explore their impact on the Tevatron 
data. Clearly, a detailed treatment of Vt) at the quantum level is also 

mandatory to perform this analysis in a consistent way. 

We remark that although CLEO data on BR{b —> sy) could preclude the existence of 
a light charged Higgs boson [0 - thus barring the possibility of the top quark decaying 
into it - this assertion is not completely general and, moreover, needs further experimental 
conhrmation. In fact, there is no direct experiment (at the level of the Tevatron anal¬ 
ysis under consideration) supporting the indirect implications on charged higgses from 
radiative B-meson decays. Originally, the bounds from CLEO data were based on the 
computation up to leading order (LO) of BR{b sy)]^. However, to this order the 


theoretical result suffered from very large uncertainties ||^ . Recently the next to leading 
order (NLO) calculation has become available 130. and the theoretical situation seems 
to be settling. The NLO calculation shows that Type I charged Higgs bosons masses are 
not restricted by 6 —>■ s y decay data either because of falling inside the experimental band 
or because of being not reliable. As for Type H charged Higgs bosons, a lower bound of 
~ 255 GeV, with an error of at most several tens of GeV, has been achieved using the 
conservative (95% C.L.) CLEO allowed band BR{b —»• sy) = (1.0 — 4.2) x 10“"^ |T^ , 


Nev¬ 


ertheless, as stated, the experimental situation is not completely settled. Recently CLEO 
has presented the preliminary new result BR{b —»• sy) = (3.15 ±0.35 ±0.32 ±0.26) x 10“^ 
which modifies the upper limit above to 4.5 x 10“^ and, therefore, it weakens the 
previous bound on the charged Higgs mass. On the other hand the ALEPH result is 


'^See and references therein. Because of the approximations used, neither of these references was 
really sensitive to the potentially large quantum effects reported here at low and high tan p. 
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BR{b —>• 57 ) = (3.11 ±0.80 ±0.72) x 10“^ [|^ which implies an upper limit (90% C.L.) of 
4.9 X 10“^. Although both results are fully compatible, the latter entails a lower bound on 
the charged Higgs mass of ~ 150 GeV at large tan/? |^-thus allowing t —> b also for 
Type II models. It is our aim to investigate, independent of and complementary to the 
indirect constraints, the decay t —> H~^ b in general 2HDM’s (Types I and II) by strictly 
taking into consideration the direct data from Tevatron on equal footing as in Ref. |^. 
This study is complementary to the supersymmetric one in Ref.^] and it should be useful 
to distinguish the kind of quantum effects expected in general 2HDM’s as compared to 
those foreseen within the context of the MSSM. 


The one-loop Feynman diagrams contributing to the decay t —>■ H~^b under consider¬ 
ation can be seen in Ref. 0: Fig. 3 (all diagrams). Fig. 4 (diagrams Cbs, Cm, Cts, ±** 4 ), 
Fig. 5 (diagram Chi) and Fig. 6 (diagram Cmi) of that reference. It goes without say¬ 
ing that the calculation of these diagrams in general 2HDM’s is different from that in 
Ref. [0, and this is so even for the Type II case since some of the Higgs boson Feynman 
rules for supersymmetric models P] cannot be borrowed without a careful adaptation of 
the couplings 0. The interaction Lagrangian describing the iff6-vertex in Type-j 2HDM 
(j = /,//) is: 

Cnlb = J,. H~ b [mt cot /? Pr ± rub aj Pl] t ± h.c. (4) 

V IMy/ 

where we have introduced the parameter aj with ai = — cot /?, an = ± tan (3. From 
the interaction Lagrangian (^) it is patent that for Type I models the branching ratios 
BR{t —> b) and BR{H^ are relevant only at low tan /?, whereas for Type H 

models the former branching ratio can be important both at low and high tan (3 and the 
latter is only signihcant at high values of tan/?. 


The renormalization procedure required for the one-loop amplitude closely follows 
that of Ref. 0 . For Type H models the one-loop counterterm and vertex structures are 
formally as in Ref. 0], whereas for Type I there are some differences. Nonetheless the 
two types of 2HDM’s can be treated simultaneously within a unified formalism as follows: 
The counterterm Lagrangian for each 2HDM model j = I, II reads 


ArW _ 
^^Hbt — 


72 M, 


H-b 


w 


SCj^^ mt cot /? Pr ± 5C^P mb aj Pl t + h.c., 


(5) 


with 


SC'i' = 


6mt 


mt 

(j) _ bmb 


'iC-6Z„.R-SZlR-6Z‘ 


Stanfl 

R- - - o -^ " ^HW f /? ) 




mb 


6v 

V 


1 

2 


1 

2 


tan/? 
6 tan /? 


- — ± - 6Zh+ + -6Z\ + - 6Z’I t r - 6Z 


tan /? 


HW 


( 6 ) 


^We have generated a fully consistent set. In part they can be found in and references therein. 
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where in the last expression the upper minus sign applies to Type I models and the lower 
plus sign to Type II ~ hereafter we will adopt this convention. 

The counterterm 5tan/?/tan/? is dehned in such a way that it absorbs the one-loop 
contribution to the decay width r(i^+ —> yielding 


6 tan/? 
tan (3 


- - Uzh± + 5Zhw- + 

V I a,' 


(7) 


The quantity 

A?) = - \SZI’ - Uzi - pfl , (8) 

niT- 2 2 

contains the (hnite) process-dependent part of the counterterm, where T/ comprises the 
complete set of one-particle-irreducible three-point functions of the charged Higgs decay 
into Vt-. Other equivalent dehnitions of this counterterm are possible, but this one, 
apart from having a clear physical meaning, automatically incorporates the corrections 
that would arise in T{H^ z/,-). 

Substituting (|^) into @ one hnally gets for the Type-j model 


sc‘i^ 


26 jii 

lsz„*-- 

+ SZhw 

tan /? ± — 


2 V 






rrit 


(9) 


We immediately see that for Type I models the one-loop correction is free of “universal” 
contributions as could be expected from our dehnition of tan /? . 


The correction to the decay width in each 2HDM can be written in the following way: 
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U) 

2HDM 


^H+b)- H+ b) 

^ H+ b) 

]\t 0 ) j\t 

^|2Re(A®)l + ^l2Re(A“)l + ;pg|2Re(AS?' + A“)] . (10) 


where the lowest-order width in the on-shell a-scheme is 


^ H+b) = 


a 


16 M2 


x 1 / 2 /t ///I ^H+ \ 

A ^ 1 , 2 ? 9/5 

^ rrit mi mi 


with 


( 11 ) 


£)( i ) _ ("^2 ^2 _ ^ q |.2 ^ _|_ ^2 ^ 2 ^ _|_ ^^2 ^2 ^ 

= {ml + ml- Mjj+)ml 
Nr = {ml + ml — Mlj+)''^'t cot^ f3 

= 2 ml ml aj cot/3 . (12) 
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The corresponding correction in the Gi^-scheme is [j^: 6{Gf) = 6 — Ar. 

The renormalized one-loop vertices for each type of model are obtained after 
adding up the counterterms (|]) to the one-loop form factors: 

Al = 6Cl + Fl 

Ar = 6Cr + Fr . (13) 

We refrain from explicitly quoting the rather lengthy expressions for the unrenormalized 
vertex functions Fr and Fr (and Ft- = FI above) for each 2HDM; the calculation and 
conventions follow those in Ref.|^. 

In the numerical analysis presented in Figs. 1-4 we have put several cuts on our set of 
inputs. From the study of the Bjorken process e'^ e~ ^ Z and the Higgs boson pair 
production e~^ e~ —> one obtains [| 

Mho + Mao >90-110 GeV , (14) 


and hence it cannot yet be excluded the possibility of a light neutral Higgs scalar, say 
below 50GeV, in general 2HDM’s. As for tan/3 we have restricted in principle to the 
segment 

0.1 < tan/3 < 60 . (15) 


For Type I models the limits are very weak while for Type H the limit at the low tan/3 
end is obtained from measurements of the process e~^ e~ ^ Z ^ /A^ 'y for t~ , light 

quarks and 6-quark decay channels . We adopt the same low tan /3 limit for Type 
I models since in this region the analysis should be similar. For the three Higgs bosons 
coupling we have imposed that they do not exceed the maximum unitarity level permitted 
for the SM three Higgs boson coupling, i.e. 


\^hhh\ ^ - ^TeV)\ - g-^ . (16) 

This condition restricts both the ranges of masses and of tan/3. Moreover, we have 
imposed that the extra induced contributions to the p parameter are bounded by the 
current experimental limit ^ : 

|Ap| < 0.003 . (17) 


With these restrictions, which are independent and truly effective in our calculation, we 
limit our numerical analysis within a wide region of parameter space where the correction 
^See Ref. for a review and references therein. 

^Notice that this condition restrains Ar within the experimental range and a fortiori the corresponding 
corrections in the Gi?-scheme. The bulk of the EW effects are contained in the non-universal corrections 
predicted in the a-scheme. 
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m itself remains perturbative, except in those places where for demonstrational purposes 
we explicitly exhibit a departure from this requirement. 


Before exploring the implications for the Tevatron analyses, we wish to show the great 
sensitivity (through quantum effects) of the decay t —>• b to the particular structure of 

the underlying 2HDM. Therefore, in the following we summarize our systematic scanning 
over the parameter space of 2HDM’s; in some cases, just to illustrate maximum effects, we 
have stretched their ranges to the very limits dehned by conditions ([T^) - ([T7|) . In all cases 
we present our results in a signihcant region of the parameter space where the branching 
ratios BR{t —>■ b) and BR{H^ —^ Ur) are expected to be sizeable. This entails 
relatively light charged Higgs bosons {Mh+ ~ 150 GeV) and a low (high) value of tan [3 
for Type I (II) models. 

In Fig. 1 we display the evolution of the correction (|l0|) with tan/? for Types I and 
II 2HDM’s and for two sets of inputs A and B for each model. We separately show the 
(leading) EW contribution, and the total correction, ^Totai = ^ew + <^qcd, which 

incorporates the conventional QCD effects 0. In this hgure we have skipped the interval 
2 < tan /? < 10 where the branching ratio of f b is too small to be of phenomeno¬ 

logical interest. In the relevant tan/? segments, that is below and above the uninteresting 
one, we hnd that the pure EW contributions can be rather large, to wit: For Type I mod¬ 
els, the positive effects can reach ~ 30%, while the negative contributions may increase 
‘arbitrarily’ - thus effectively enhancing to a great extent the modest QCD corrections- 
still in a region of parameter space respecting the restrictions 0 - 0 ; For Type II mod¬ 
els, instead, the EW effects can be very large, for both signs, in the high tan /? regime. 
In particular, the huge positive yields could go into a complete “screening” of the QCD 
corrections. 

In Fig. 2 (resp. 3) we present the evolution of the corrections in Type I (resp. II) models 
as a function of the other parameters, namely the tangent of the CP-even mixing angle 
(a), the charged Higgs mass (b), the CP-odd scalar mass (c) and the CP-even scalar 
masses (d). Inputs A and B for each model are as in Fig. 1 whenever they are hxed. 
Notice the rapid oscillation, yet qualitatively different for both models, around tana = 0. 
On the other hand, the fast evolution that can be seen as a function of the masses is due 
to the fact that the Higgs self-couplings in the three-point functions are proportional to 
the splitting of the Higgs masses. In Fig. 2d (also in Fig. 3d) the range of the CP-even 
masses is plotted until the condition M^o < M^o is exhausted or there is a breakdown of 
relations ([IBD and/or (0)- 


Next we turn to the discussion of the dramatic implications that the EW effects may 
have for the decay t b at the Tevatron. The original analysis of the data (based 

on the non-observation of any excess of r-events) and its interpretation in terms of limits 
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on the 2HDM parameter space was performed in Ref. IQ (for Type II models) withont 
inclnding the EW corrections. In these references an exclnsion plot is presented in the 
(tan/?, Mji^+)-plane after correcting for QCD effects only. To demonstrate the potential 
impact of the EW loops on these studies we follow the method of Ref. 0. Although the 
data used by the Tevatron collaborations is based on inclusive r-lepton tagging |^, it 
will suffice for illustrative purposes to concentrate ourselves on the (r,//-channel ||, . 


In this way the comparison of the results for generic Type II 2HDM’s and those already 
available for the specihc case of the MSSM Higgs sector |p will be more transparent. The 
production cross-section of the top quark in the (r,//-channel can be easily related to the 
decay rate of / ^ b and the branching ratio of Uj. as follows: 


— 


81 ^ 9 r(/ 


4 r(/ ^ H+b) 
IT+ b) 


BR{H+ T+ Vr) €2 




(18) 


This formula generalizes eq.(7) of Ref. for the case that BR{H —> TUr) is not 100%, 
as it indeed happens when we explore the low tan (3 region. In general we have 

T{H+ T+ Ur) 


BR{H^ 


Ur) = 


T{H+ r+ Ur) + T{H+ ^ cs) ’ 


(19) 


where we use the QCD-corrected amplitude for the last term in the denominator . 
This branching ratio is about 50% for Type I 2HDM at low tan /?, and 100% for Type II 
at high tan/? (the case studied in P). 


Finally in Figs. 4a and 4b we have plotted the perturbative exclusion regions in the 
parameter space (tan/?, for intermediate and extreme sets of 2HDM inputs A, B, 

B’ and C. In Type I models (a) we see that the bounds obtained from the EW-corrected 
amplitude are generally less restrictive than those obtained by means of tree-level and 
QCD-corrected amplitudes. Evolution of the excluded region from set A to set C in 
Fig. 4a shows that the region tends to evanesce, which is indeed the case when we further 
increase M^o in set C. In Type II models (b) we also show a series of possible scenarios. We 
have checked that the maximum positive effect 6ew > 0 (set A in Fig.4b) may completely 
cancel the QCD corrections and restore the full one-loop width —> H~^ b) to the tree- 

level value ( pTl) just as if there were no QCD corrections at all! Intermediate possibilities 
(set B’) are also shown. In the other extreme the (negative) effects 6ew < 0 enforce 
the exclusion region to draw back to curve C where it starts to gradually disappear into 
a non-perturbative corner of the parameter space where one cannot claim any bound 
whatsoever!!. 


Some discussion may be necessary to compare the present analysis with the super- 
symmetric one in Refs. 0, |^. In the MSSM case, the Higgs sector is of Type H. However, 
due to supersymmetric restrictions in the structure of the Higgs potential, there are large 








cancellations between the one-particle-irreducible vertex functions, so that the overall 
contribution from the MSSM Higgs sector to the correction (|I0|) is negligible. In fact, we 
have checked that when we take the Higgs boson masses as they are correlated by the 
MSSM we recover the previous result 0. Still, in the SUSY case there emerges a large 
effect from the genuine sparticle sector, mainly from the SUSY-QCD contributions to the 
bottom mass renormalization counterterm [^] , which can be positive or negative because 
the correction flips sign with the higgsino mixing parameter. In contrast, for general 
(non-SUSY) Type H models the bulk of the EW correction comes from large unbalanced 
contributions from the vertex functions, which can also flip sign with tana (Cf. Fig. 3a) - 
a free parameter in the non-supersymmetric case. Although the size and sign of the effects 
can be similar for a general Type H and a SUSY 2HDM, they should be distinguishable 
since the large corrections are attained for very different values of the Higgs boson masses. 
For instance, in generic 2HDM’s (of both types) large negative effects may occur for large 
values of the CP-odd Higgs mass (Cf. Figs. 2c and 3c). In the MSSM the latter should 
be essentially degenerate with the charged Higgs in that region and so t b would 

have never occurred. 

Therefore, just to illustrate one possibility, let us envision the following scenario. Sup¬ 
pose that t —>■ b is not observed at the Tevatron - or that it comes out highly sup¬ 
pressed beyond QCD expectations (Cf. Fig.4b, curve C) - while at the same time H~^ 
and are observed (maybe produced at the Tevatron itself or at LEP) and both show 
up with a similar mass below m*. Then these bosons could well be supersymmetric Higgs 
bosons. If, however, a similar situation would be encountered but A° is not produced (be¬ 
cause it is perhaps too heavy), then the observed cannot probably be a SUSY Higgs. 
Notice (Cf. Fig. 4a) that in this case the could also belong to a Type I model. In this 
case further investigation would be required to disentangle the type of non-SUSY Higgs 
model at hand e.g. aiming at a determination of tan/?. One possibility would be from 
the decay H'^ —> after been produced from mechanisms other than top quark 

decay; alternatively, once a heavy A° would be found it would provide a handle to a tan (3 
measurement through the decays — >■ t~ and/or A^ ^ bb. If these decays would be 
tagged at a high rate, the non-SUSY model should necessarily be of Type H. At the LHC, 
or at a NLC, one could even use P(A° —> bb)/r{A^ tt) (x tan^/? for very heavy A^. In 
short, a combined procedure based on quantum effects and direct production could be a 
suitable strategy to unravel the identity of the Higgs bosons. 

To summarize, we have computed the electroweak one-loop corrections to the uncon¬ 
ventional top quark decay width P(f —> b) at the leading order in the Yukawa couplings 

both for general Type I and Type H 2HDM’s. We have found that the EW corrections 
can be comparable in size to the QCD effects and be of both signs. The positive ones can 
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reach 30% and 50% for Type I and Type II models respectively, which means that they 
conld simply delete the QCD corrections. The negative ones can even be larger (—50% 
or more in ample regions of parameter space) for both models. We have also shown that 
this fact may deeply inflnence the cnrrent interpretation of the Tevatron data on that 
decay. Most important, we have argned that knowledge of the EW qnantnm effects may 
be crncial to nnderstand the natnre (snpersymmetric or else) of the Higgs bosons, if they 
are eventnally fonnd in fntnre experiments at hadron and/or e"*" e“ colliders. 
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Figure Captions 

• Fig. 1 The correction 6, eq. (0. to the decay width P(t —> LL~^b) as a function 
of tan/3, for Type I 2HDM’s (left hand side of the figure) and two sets of inputs 
{{Mh+, Mfjo, Mho, Mao); tana}, namely set A: {(70,175,100, 50) GeV; 3} and set 
B: {(120, 200, 80, 250) GeV; 1}. Similarly for Type II models (right hand side of the 
hgure) and for two different sets of inputs, set A: {(120, 300, 50, 225) GeV; 1} and 
set B: {(120, 300, 80, 225) GeV; —3}. Shown are the electroweak contribution Jew 
and the total correction Jxotai = Jew + Jqcd- 

• Fig. 2 The corrections Jew and Jxotai for the Type I 2HDM as a function of (a) 
tana, (b) the charged Higgs mass, (c) the pseudoscalar Higgs mass, and (d) the 
heavy (labelled with a triangle) and light (unlabelled) scalar Higgs masses. Inputs 
as in Fig. 1 with tan/3 = 0.1 for set A and tan/3 = 0.2 for set B. 

• Fig. 3 The corrections Jew and Jxotai for the Type H 2HDM as a function of (a) 
tana, (b) the charged Higgs mass, (c) the pseudoscalar Higgs mass, and (d) the 
heavy (labelled with a triangle) and light (unlabelled) scalar Higgs masses. Inputs 
as in Fig. 1 with tan/3 = 35 for both sets. 

• Fig. 4 The 95% G.L. exclusion plot in the (tan/3, Mji/+)-pIane for (a) Type I 2HDM 
using three sets of inputs: A and B as in Fig. 1, and G: {{Mh+, 200, 80, 700) GeV; 1}; 
(b) Similarly for Type H models including three sets of inputs: A as defined in Fig. 1, 
B’:{(Mjy+, 200, 80,150) GeV; 0.3}, and G:{(Mh+, 200, 80,150) GeV; -3}. Shown 
are the tree-level, QGD-corrected and fully 2HDM-corrected contour lines. The 
excluded region in each case is the one lying below these curves. 
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